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Abstract 

We study the exclusive decays of B — ► fo(980)K^ in the framework of the perturbative QCD 
by identifying the /o(980) as the composition of ss and nn = (uu + dd)/y/2. We find that the 
influence of the nn content on the predicted branching ratios is crucial. We discuss the possible 
rescattering and gluonium states which could enhance the branching ratios of considered decays. 
We point out that the CP asymmetry in B — > fo(980)Ks,L could be a new explorer of sin2</>i. 
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Recently, Belle |l| and BaBar 0] have opened new channels in three-body nonlep- 
tonic B decays, such as B — > KKK, KKn and Knir decays. In particular, BELLE 
has observed the decay of B + — > / (980)-ft' + with the branching ratio (BR) product of 
Br(B+ - /o(980)ir+) x 5r(/o(980) -> tt+tt-) = (9.6±i : |±i;i±g:|) x lO" 6 0. Since / (980) is 
a neutral scalar meson, the measured B — > /o(980)i^ decays not only show that in the first 
time -B decays to scalar-pseudoscalar final states but also provide the chance to understand 
the content of / (980) and its production. 

The essential inner structure of the scalar meson / (980) is still obscure since it was 
established first by Ref. [[§ with phase shift analysis. In the literature, /o(980) could be 
four-quark states denoted by qqqq [§J or KK molecular states || or qq states . However, 
one objection against the possibility of KK states is that the binding energy of 10 — 20 MeV 
for KK is much smaller than measured width in the range 40 — 100 MeV 0. It is suggested 
that in terms of the measured — > f (980)7 an d /o(980) — > 77 |7|, @, H [HJ decays and 
— > / (980)7T + decay ]10|, |ll|], the flavor contents of /o(980) are ss mostly and a small 
portion of nn = (uu + dd)/y/2. In this paper, we take /o(980) to be composed of qq states 
and use |/o(980) >= cos0 s |ss > +sin0 s |fm > to denote its flavor wave function. 

Before presenting our perturbative QCD (PQCD) calculation to the decays, we would 
like to give a brief model-independent analysis on B — > fo(980)K. For simplicity, our 
analysis will only concentrate on dominant factorizable parts and regard the /o(980) as the 
composition of ss so that at the quark level the process corresponds to b — > sss decay. Since 
at the quark level, the physics for B — > <j)K decays is the same as B — ► /o(980)ii^*), we 
first examine the <j)K mode. We start by writing the effective Hamiltonian for the b — *■ s 
transition as 
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q'= u ,, 



i=3 

where V g / = V*, s V q >b are the products of the CKM matrix elements, Cj(/i) are the Wilson 
coefficients (WCs) and Oi correspond to the four-quark operators. The explicit expressions 
of Ci(fi) and Oi can be found in Ref. \Y1 . It is known that the vector meson has the 



current matrix elements 

(O|s7 M s|0) = m^e^, (0\ss\<f>) = 0. (1) 
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Therefore, the decay amplitude of B d — > 4>K° can be simply described by JT3 

A (B d -> <f>K°) = UV t * (af + a[ s) + a { 5 s) ) F e BK + 2f B V t *af F * 6 K + ... (2) 
where af^ are defined by 

n i ^2 Ci 
a 3,4 — ^3,4 + —^9,10 + 0-3,4 5 a 3,4 ~ ~J^~ + ^-^10,9, 



l 5,6 - ^5,6 + ~~^~ U 7,8 + «5,6 > a 5,6 ~ + 2N~ J ' 



Note that a 2 is larger than a\ and a\ 6 are much larger than af\, Ff K is the B —> K decay 
form factor and F^f denotes {(^K^^d^S) annihilation effect. We neglect {(f)K\s^ ^<jd\Q) 



due to the chirality suppression. According to the results of Ref. [JL4J] , the predicted BR of 
B d — > <pK° is around ~ 9.6 x 10~ 6 which is consistent with Belle's and BaBar's results given 



by (10.0±1;?±?;§) x 10~ 6 U and (8.7i{;£ ± 0.9) x 10~ 6 0, respectively. However, in the 



7+1.7 



B — > f (980)K° decay the relevant current matrix elements are given by 

(0|*y„s|/ (980)) = 0; (0|ss|/ (980)> = m f J (4) 

where rrif (f) are the mass ( decay constant) of /o(980). From Eqs. ([I]) and (|), it is 
clear that the situation in the f (980)K° mode is just opposite to the case of 4>K°, i.e., the 
role in vector and scalar vertex is exchanged each other. Hence, the decay amplitude for 
B — > f (980)K° can be simply written as 

A (B d - / (980)K°) = 2r f f~vM s) Sf K + 2f B V?a!f>F£ K + ... (5) 

where r/ = mf /MB, = (K\bs\B) /Mb and the factor of 2 comes from Fierz transforma- 
tion in the (V — A) x (V+A) operator. By taking Mb ~ raj, with mj being the 6-quark mass, 
one can expect S^ K ~ Ff K by equation of motion. Comparing Eq. (§) with Eq. (H), it is 
obvious that there is a suppression factor 77 = 0.186 in B — > /o^SO)^ . Including the fac- 



tor of 2, taking a§\\J KMb) / af\\j KMb) ~ 1.5 and / ~ and neglecting the annihilation 
contributions, we estimate the ratio of BriB — > fo(980) K°) / Br (B — > 0if°) being around 
0.31. With the average value of Belle and BarBar, the predicted BriB — ► /o^SO)^ ) is 
around 2.8 x 10~ 6 . Will the value change while we include the nn content and annihilation 
contributions? In the following we display the results in a more serious theoretical approach. 



It is known that large uncertainties are always involved in the calculations of transition 
matrix elements while studying exclusive hadron decays. Nevertheless, the problem will 
become mildly in the heavy B meson decays because of the enormous B samples produced 
by B factories. That is, more precise measurements will help the theory to pin down 
the unknown parameters to make some predictions. In this paper, we adopt the PQCD 
approach in which the applications to exclusive heavy B meson decays, such as B — > |T9| , 
B -> mr(KK) |(J |2|, B -> <fm(K) 0, g|], B -> rj^K H and B -+ pK p| decays, 
have been studied and found that all of them are consistent with the current experimental 



data p5| , ^q] . In the PQCD, in order to solve the various divergences encountered at end- 
point, we will include not only resummation, for removing end-point singularities, but 
also threshold resummation, for smearing the double logarithmic divergence arisen from 
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In order to satisfy the local current matrix elements with Eq. (Q), the /o(980) meson 
distribution amplitude is given by 

(0|g(0) iff (z)»|/ o > = -^J* dxe-^mfilj^f^x) (6) 



with the normalization 



dx® fo (x) 



lo Joy ' 2^2N~ C 

For B and mesons, the corresponding distribution amplitudes |14|, [17| are written as 

1 f 1 



(0|6(z) ig (0)i|S> = -/=y o dxe~ inP ''{[JP] jl + M B [r\ jl )>y B Mz) (7) 
(Jjr|g(z) ia (0),|0> = ^^W^{[ 75 f\j&K{x) 

+rn° K hs} j i$ p K (x) + [75 (A- " (8) 

(!T( eji )|gWia(0),|0> = / Q W- p -{Mk4M;^*(*) 

} (9) 

with n_ = (0, 1, 0j_), n + = (1, 0, 0j_) and m° K being the chiral symmetry breaking parameter, 
$ x m (x) stand for twist-2 wave functions and the remains belong to twist-3 with their explicit 



expressions being given in Ref. ||14| , |lq| . Although K* meson has three polarizations, only 
the longitudinal part is involved in the decays of B — > /o (980) if*. We note that as the K 
meson, the distribution amplitude of / (980) can have more complicated spin structures. 



However, since our purpose is just for the properties of B — > f (980)K ( -*\ we consider only 
the simplest case. Moreover, if / (980) consists of ss mostly, the choice of Eq. @ is clearly 
dominant . 

The decay rates of B — > f (98Q)K are expressed by 



G 2 F M s BlAl2 



32tt 



A 



(10) 



where A includes all possible components of /o(980) and topologies. As mentioned before, 
/o(980) has the components of ss and nn, for different contents, the amplitudes of B d — > 
/ (980)i?° and B + -> / (980) K + are written as 



Az 



Af. 



A$s 
A+ 



f~v t *s2 s) + f B v:s p Ji ] + 



f K V t *N^ d) + f B V t *N^ + ..., 

fVt^ds^ + IbV^S^ - f B V*S a + 

f K V t *N^ + f B V t *N^ ] - f K V:N e - f B V:N a + 



(11) 



respectively, where S^ffl (^e(a)) denote the emission (annihilation) contributions of the 
ss (nn) content from penguin diagrams while S e ( a )(N e t a )) are from tree contributions. The 
total decay amplitude for the neutral (charged) mode is described by A = cos 4>A^ + 
sin 4>A^n /V2- Due to the smallness of nonfactorizable effects, we neglect to show them in 
the Eq. ( pTTj) and just display the factorizable contributions with emission and annihilation 
topologies. But we will include their effects in our final numerical results. According to Eqs. 
(0— Bl), the hard amplitudes for ss content are derived as 



S. 



e6 



16nr (Ol s CfM\ \ dx\dx 3 \ b\db\b 3 db 3 ^ B (x\,b\) 
Jo Jo 

x{[^ K {x 3 ) + 2r K ^ p K {x 3 )+r K x 3 {^ p K {x 3 ) - $£(x 3 )) 

xES^{t^)h e (x 1)X3 MM) 

+2r K ^ K (x 3 )ESif(t^)h e (x 3 ,x h b 3 ,b 1 )}, 



(12) 



°a4 



rl roo 

-lQnr fr K C F M B \ dx 2 dx 3 / b 2 db 2 b 3 db 3 § r (x 2 ) 
Jo Jo 

"(1 + x 3 )& K (l - x 3 ) + (1 - x 3 W K {\ - x 3 )] E$(tW)h a (x 2 , x 3 , h, h) 
l + x 2 )^ p K (l-x 3 )}E^(t^)h a (x 3 ,x 2 ,b 3 ,b 2 )\, (13) 
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S^ q) = 16nr f C F M 2 B J dx 2 dx 3 J™ b 2 db 2 b 3 db 3 <f> fo (x 2 )$ K (l - x 3 ) 

x \2E { ^(t^)h a (x 2 , x 3 , b 2 , b 3 ) + x 2 E ( £(t^)h a (x 3 , x 2 , 63, b 2 )\ (14) 



and the results for the content of nn are given by 

N^ q) = 8nr f a s C F M B J dx x dx 2 J b 1 db 1 b 2 db 2 §B(xi, bi)$f (x 2 ) 

x {(1 - 2x 2 )EN^(t^)h e (x 1 , x 2 , b u b 2 ) + 2ENi}(t^)h e (x 2 , x u b 2 , h)} (15) 

NeQ 9 ^ = —I6nr frKCXsCpMg J dx\dx 2 J bidbib 2 db 2 & B (x\, bi)&f (x 2 ), 

x {(2 + x^EN^^h^, x 2 , b u b 2 ) + 2EN^(t^)h e (x 2 , x u b 2 , 6 X )}. (16) 

Here, the hard part functions h e ^ , mainly arising from the propagators of hard gluon and 
valence quark, have included the threshold resummation factor [|rjj] , and the evolution factors 
are given by 

ES { £(t) = at a (t)a$\t)exp[-SB(t,x 1 )-S K (t,x 3 )], 
EN { j\t) = aL a (t)aW(t)exp[-S B (t,x 1 )-Sf a (t,x 2 )], 
E%\t) = a s (t)a^\t)ex P [-S fQ (t,x 2 ) - S K (t,x 3 )\ 

where the exponential effects denote Sudakov factors generated by the k F resummation. 
Since the annihilation topologies in both contents are associated with the matrix elements 
< f (980)K\s>y u >y 5 d(u)\0 > and < f (980)K\s'y 5 d(u)\0 >, the difference in the different 
content is only from the spectator, which is s (d or u) in ss (nn) component. Hence, 
could be obtained from S^l by replacing x 2 with 1 — x 2 in &f (x 2 ) and 1 — x 3 with x 3 in 
{$ K (1 - x 3 )}. In Eq. (0), we define that = S^ q) + S^ q) and = N*} q) + N^ {q) . 

The S a and N e can be obtained from S^} 9 ^ and N^ 9 ^ by replacing W.C. with a 2 . We note 
that unlike the cases of B — > PP and B — > PV decays in which the chirality suppression in 
the (V — A) x (V — A) annihilation topology is conspicuous, the contribution of Eq. ([13]) 
may not be small because one twist-3 effect, such as is not cancelled. And also, differing 
from B — > PP decays, Eq. (|l|) is opposite to Eq. (|T6| ) in sign. 

Similar to B -> / (980) if decays, the amplitudes for B -»■ f (980) K*° and B + -»■ 
/o(980)K* + can be expressed as 

A- ss = fV t *K2 s) + f B V t *KS 6 d) + ..., 

Ann = fK*V t *X^ d + f B V f *X^Q + 



At 



f~V t *KS s) + f B V t *K^ - f B V*K a + 

f K *v;xS u) + f B v;x^ - f K *v:x e - f B v*x a + .. 



(17) 



where K^ q \ Ka} and can be obtained from S^Q q \ and S^ q \ respectively, by 

replacing (fix, 4> P k an d <Pk with 4>k*, 4>k* an d <$k*- With the similar analysis, we find that 

X e 4 = 



iV e 4 and X^l are related to K^jL if we change £2 and 1 — £3 to 1 — £2 and £3, 



L a4(6) 



respectively. The -K" a and X e are the same as K a ^ and X e4 but the associated WC is a 2 . 



In our numerical calculations, the B meson wave function is taken as [2~I| 



$ B (x,b) = N B x 2 (l - x) 2 exp 



where uj b is the shape parameter and N B is the normalization, determined by 

Sb 





^ xM B N 


) 2 


-5 







/ dx$ B (x,6 = 0) 

JO 



Since /o(980) is a light meson, its wave function can be defined in the frame of the light-cone 
and the concept of the twist expansion can be used. And because the relevant scalar meson 
wave function has not been derived in the literature yet, we choose the following form 

$ /0 (£) = ^^{3(l-2x) 2 + e(l-2£) 2 (C' 2 3/2 (l-2£)-3) + 1.8C' 4 1/2 (l-2x)} 

with Cl /2 (y) = (35y 4 - 30y 2 + 3)/8, Cl /2 (y) = 3/2(5y 2 - 1) and f = 0.3 ~ 0.5 in our 
estimations. In order to fix the values of u B and m° K , we take uo B = 0.4 and m° K = 1.7 



GeV as those in the studies of the B — ► K and B —> K* form factors [21, 27]. Explicitly, 
with the values above and / = 0.2 0, the B — > /o(980) form factor is found to be 0.270 
(0.286) for £ = 0.3 (0.5). Using Eqs. (11 — 17), the values of hard amplitudes are shown in 
Table |. In addition, the decay BRs and CP asymmetries (Acp) for B — > / (980)if and 
B -> /o (980) if* are shown in Table |TJ with <p s = and in Table [TTI| with <p s = 40° and 140°, 
where 

r-r 



A 



CP 



r + r 



with T (f ) being the particle (antiparticle) partial decay rate. In both tables, the CP 
violating phase of V u b is taken to be <p 3 = 72°. From Table |, we see clearly that the sign 
of annihilation contribution S^fg is opposite to so that the BRs of B — > f (980)K 

are reduced significantly. Because the corresponding value of S^}^ in B — > f (980)K* is 



TABLE I: The parameters in the amplitudes with ujb = 0.4, m\ = 1.7 and / = 0.2 GeV for 
i = 0.3. 



Amp. 


s£ w (io- 2 ) 


sffi (io- 2 ) 


5„(10- 2 ) 


iv2 u) (io- 2 ) 


iVe(10- 2 ) 


/ (980)K 


-1.02 


0.45 + i0.39 


-7.41 - i0.40 


1.75 


31.44 


Amp. 


< 6 (s) (io- 2 ) 


<4 ( 6 d) (10- 2 ) 


^a(10" 2 ) 


X e P 4 (W) (lO- 2 ) 


x e (io- 2 ) 


/ (980)K* 


-1.37 


0.06 + i0.40 


0.83 - il.98 


1.31 


-31.44 



too small as shown in Table ||, we can understand that the influence of the annihilation 
is insignificant. Also we see that if / (980) is ss mostly, BRs of B — > f (980)K are only 
around 10 -6 . On the contrary, they could be over 5.0 x 10~ 6 at a large angle of <f) s ; whereas 
BRs for B — > /o(980)i^* are similar to those at <p s = 0. To further understand the effect of 
nn, we display the BRs as a function of <p s in Figure [I]. Obviously, the predictions are very 
sensitive to the contributions of the nn component. The essential question is what the rang 



of <p s is allowed. Unfortunately, the preferable 4> s is still unknown, i.e., both small |28[] and 
large |§ angle solutions exist simultaneously in the literature. The former prefers the value 
of S = 42.14i|;|° but the latter 138° ± 6°. 

It is also interesting to look for a decay which is directly related to the nn content only 
such that it can tell us how the impact of the content would be. We point out that one 
of the possible decays is B + — > / (980)7r + . In this decay, from the identity of Eq. (|J) 
the contribution from the ss content vanishes since it corresponds to the vector current 
vertex. As the 7r meson is similar to K except 577(3) breaking effects as well as the relevant 
wave functions and CKM matrix elements, the BR is estimated to be 0.25 (0.45) x 10~ 6 for 
(f) s = 35° (132°). 

It seems that in spite of the uncertain allowed value of S , the predicted BR of B + — > 
/o(980)K + is smaller than the central (minimal) value of 21.17 (11.03) x 10~ 6 reported by 
Belle in which 5r(/ (980) -> tt+tt-) = 2/3R with R = r(7T7r)/r(7r7r) + T(KK) ~ 0.68 |2j 
is used. It is clear that unless there exists some other mechanism in the decay processes, it 
is difficult to explain the large measured BR by just considering short-distant interactions. 
Actually, the case is similar to the charmed decay of B + — > XcoK + with its BR being 
(6.0i 2 ;g±l.l) x 10~ 4 m and (2.4±0.7) x 10~ 4 measured by Belle and BaBar, respectively. 
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FIG. 1: BRs of (a) B -» /o(980)A and (b) B -» / (980)ifM as a function of S . The bold (thin) 
solid lines stand for neutral modes with £ = 0.3(0.5) while the dashed lines are for charged modes. 

TABLE II: BRs of B — » /o(980)i^M without and with annihilation contributions by taking the 
same values of parameters as Table |. The CPAs with annihilation effects are also shown. 





Mode 


BR(xl0~ 6 ) 
no anni. 


BR(xlO~ 6 ) 


CPA (%) 


B - 


-> /o(980)^ 


2.95 


1.39 





B+ - 


- / (980)A'+ 


3.13 


1.57 


6.50 


B - 


* /o(980)A*° 


5.16 


5.40 





B+ - 


- /o(980)A*+ 


5.49 


5.76 


1.48 



In order to to agree with the experimental data, the authors of Ref. |32|] suggest that the 
possible rescattering effects, such as B — »■ D*D, B — > D S D* and B — > D*D* decays via 
triangle diagrams with the strong couplings of <?dd Xc0 , gD*D* X co> 9d*d* X co e ^ c -; could enhance 
the BR and reach (1.1 ~ 3.2) x 10~ 4 . We find that the mechanism could also apply to 
B — > fo(980)K^ decays by replacing Xco with / (980) and taking the proper couplings such 
as 9D*D s f (980) etc.. The illustrated diagrams are displayed in Figure |2|. Since the differences 
between f (980)K^*> and XcqK + modes are only related to the strong coupling constants, by 
following the procedure of Ref. |J2] and taking proper values of strong coupling constants, 
one expects that the rescattering effects in Figure can reach (9(10 -5 ) easily. In addition, 
we note that, as the possibility of gluonium state in rj' W\ B3J was proposed to explain the 
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TABLE III: BRs and CPAs of B -> / (980)if W decays with the same parameters as Table § but 
by taking cp s = 40° and 140° . 



Mode 




BR,(xl(r 6 ) 


CPA (%) 


B - 


-> /o(980)^° 


40° 


1.14 









140° 


4.70 





B+ - 


-> /o(980)K+ 


40° 


1.25 


-16.94 






140° 


5.94 


1.58 


B - 


♦ /o(980)iT Q 


40° 


0.85 









140° 


6.70 





B+ - 


-» /o(980)if*+ 


40° 


1.51 


11.31 






140° 


6.30 


-12.19 




FIG. 2: Possible rescattering diagrams contributing to B -> f (980)K^ 



tendency of large BR for rf production fl34~l , it is also possible to understand B — > / (980).K" 



decays by considering g — g — /o(980) coupling. 

Finally, it is worth to mention that like the decays of B — > J/^Ks^l and S — > </)Ks,l, 
B — > /o(980)i ; fs i L can be another outstanding candidate for the observation of the time- 
dependent CP asymmetry, defined by 

F(t)-r(t) 



r(*) + r(t) 

with 



-> /) cos(AM d t) - A^(B d - /) sin(AM d ) (18) 



1 + |A. 



2 
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where / expresses the CP final state classified by n = ±1, A^'p (A™™) denotes the direct 
(mixing- induced) CP violation and A(B(B^) — > /) are the B(B) meson decay amplitudes 
and 0i = Arg(V t * d V t b)- From Eq. (0), we know that if no CP violating phase is involved 
in the decay amplitude, the direct CP violating observable will vanish and A = rje' 2 " 1 ^ 1 . 
Consequently, the mixing- induced CP asymmetry is only related to sin20i. Fortunately, 
since there is no tree contribution, B — > fo(980)Ks,L decays satisfy the criterion so that 
the mixing- induced CP asymmetry of B — > f (980)Ks,L is expected to be the same as that 
measured in B — ► J/^Ks } l decays, except a small derivation from higher order contributions 
3|. 



We have studied /o(980) scalar meson production in B meson decays by assuming 
that its flavor contents are ss and nn states. We have found the role of nn on the 
BRs of B — > fo(980)K is crucial. We have also pointed out that nonvanishing BR of 
B + — ► / (980)7T + could test the existence of the nn content although the expected BR is 
less than 10" 6 . If the BRs of B -> f (980)K^ are measured to be 10 5 in experiments, 
the results could be the evidence of existing rescattering effects in heavy B meson decays. 
On the other hand, although we concentrate on the qq states of / (980), we do not 
exclude the possibilities of four-quark qqqq and gluonium states. The possibilities could 
be clarified by other experiments, such as — > Tm^KK)^ D+ — > 7r + 7r + 7r~ |J7| decays 
and B + — > / (980)7r + with the measured BR of O(10~ 6 ). Furthermore, the more general 



approach to deal with the decays of B — ► f (980)K can refer to Ref. ||38| . We have also 
suggested that time- dependent CP asymmetry in B — > 70(980)^5^ could be a new explorer 
of sin 20i. We note that our study of B decays with /o(980) in the final state can be 
applied to other modes with scalar or isoscalar final states such as a(600), ao(980) etc.. 
Remarkably, B factories have opened a new opportunity to understand the contents of 
scalars and their productions. 
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